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SECTION  I 
INTRODUCTION 


Certain  experiments  and  tests  require  the  measurement  of 
out-of-plane  motions  of  target  surfaces.  One  way  these  measure¬ 
ments  can  be  made  is  with  use  of  special  devices  which  produce 
moire  fringe  patterns  on  the  surface  of  the  targets.  A  descrip¬ 
tion  of  one  of  these  devices  has  been  provided  in  a  previous 
document.  This  document  provides  detailed  descriptions  of 
procedures  related  to  the  use  of  that  device  and  the  analysis 
of  the  resulting  fringe  patterns. 

Although  the  device  has  been  used  primarily  to  determine 

the  response  of  impulsively  loaded  targets  of  fairly  large 

2 

area  (1500  to  2000  cm  ) ,  information  provided  in  this  document 
also  applies  to  analysis  of  data  from  surfaces  of  smaller  area 
and  surfaces  subjected  to  slowly  applied  or  static  loads.  In 
all  cases,  deflection  data  are  obtained  after  analysis  of 
photographs  of  individual  fringe  patterns. 

A  computer  program  named  PANDIS  is  used  to  assist  in  the 
reduction  of  fringe  pattern  data  to  displacement  data.  This 
program,  written  in  the  BASIC  language,  is  designed  for  use  from 
a  remote  time-shared  computer  terminal.  Section  II  of  this 
document  provides  a  brief  description  and  instructions  for  use 
of  the  device.  Sections  III,  IV,  and  V  provide  detailed  data 
reduction  instructions,  a  description  of  program  PANDIS  (in¬ 
cluding  a  listing  of  program  statements) ,  and  a  sample  problem 
and  the  resulting  printout,  respectively. 


(1)  Piekutowski,  A.  J.,  A  Device  to  Determine  the  Out-of-Plane 
Displacement  of  a  Surface  Using  a  Moire  Fringe  Technique, 
AFWAL-TR-81-300 5 ,  Air  Force  Wright  Aeronautical  Laboratories, 
Wright-Patterson  Air  Force  Base,  OH,  May  1981. 
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SECTION  II 
USE  OF  MOIRE  DEVICE 

As  constructed,  the  moire  device,  shown  schematically  in 
Figure  1,  simply  generates  a  system  of  light  and  dark  bands 
which  are  fixed  in  space  relative  to  the  moire  device.  The 
light  and  dark  bands  are  apparent  to  the  viewer  only  when  the 
system  of  bands  is  interrupted  by  a  target  surface.  A  critical 
orientation  of  the  target  surface  with  respect  to  the  moire 
device  is  not  required  since  the  initial  shape  of  the  target 
surface  can  always  be  reconstructed  after  analysis  of  the  fringe 
patterns . 

An  illustration  of  the  technique  used  to  produce  the  moire 
patterns  is  presented  in  Figure  2.  As  shown  in  this  figure, 
the  optical  system  on  the  right  (Projection  System)  projects 
an  image  of  the  Ronchi  ruling  in  the  system  onto  the  target 
surface.  The  lens  in  the  left  optical  system  (Viewing  System) 
focuses  an  image  of  the  target  and  projected  ruling  onto  the 
Ronchi  ruling  in  the  viewing  system.  A  moire  pattern  is  pro¬ 
duced  where  the  image  of  the  projected  ruling  is  superimposed 
on  the  viewing  system  Ronchi  ruling.  A  camera,  with  its  lens 
focused  on  the  composite  image,  is  used  to  photographically 
record  the  moire  pattern. 

Analysis  of  raw  data  is  facilitated  by  exercising 
reasonable  care  during  setup  of  the  moire  device.  Attention 
is  drawn  to  the  relationship  of  the  shaded  regions  adjacent 
to  the  rays  (boundaries  between  clear  and  opaque  areas  of 
Ronchi  ruling)  shown  in  Figure  2.  The  regions  shown  are 
adjacent  to  the  optical  axes  of  the  projection  and  viewing 
systems.  When  the  rays  associated  with  these  regions  and  the 
target  surface  intersect  simultaneously,  a  dark  band  will 
appear  at  the  target  surface.  This  particular  dark  band  is 
identified  as  fringe  zero.  Identification  of  the  rays  which 
produce  fringe  zero  is  facilitated  by  heavy  locator  lines 


BOTH  SYSTEMS  J2Q  NOT 


NEED  TO  BE  AT  SAME  ANGLE 


TO  THE  TARGET  SURFACE 


placed  on  the  Ronchi  rulings.  Relative  positions  of  the  locator 
lines  in  the  various  dark  bands  of  the  moire  pattern  are  used  to 
establish  the  direction  which  the  moire  device  must  be  moved  to 
produce  the  desired  target-device  relationship.  The  positions 
of  the  locator  lines  are  shown  for  a  number  of  fringe  surfaces 
in  Figure  3.  The  coordinate  system  and  fringe  numbering  con¬ 
vention  which  was  adopted  for  use  during  subsequent  data  reduc¬ 
tion  procedures  are  also  shown  in  this  figure. 


PROJECTION  OF  OPTICAL  AXIS 


Figure  3.  Arbitrary  System  of  Identifying  Fringes  and  Coordi¬ 
nate  System  Used  in  Analysis  of  Fringe  Patterns.  In 
this  figure,  the  reference  plane  is  parallel  to  the 
X-D  plane.  Note  the  position  of  locator  lines  in 
each  fringe. 


Information  presented  in  Figure  4  is  used  to  illustrate 
relationships  between  the  various  device  parameters  which  must 
be  determined  for  correct  interpretation  of  the  resultant 
fringe  patterns.  However,  one  important  difference  between 
the  device  and  the  illustration  in  Figure  4  exists.  For  pur¬ 
poses  of  simplifying  presentation,  the  Ronchi  ruling  has  been 
shown  between  the  appropriate  lens  and  the  target  rather  than 
behind  the  lens,  e.g.,  Figures  1  and  2.  Positioning  the  Ronchi 
ruling  as  shown  in  Figure  4  merely  facilitates  a  graphical 

presentation  of  the  rays  used  to  determine  the  location  of  the 

dark  bands.  For  consistency  of  illustration,  the  lens  on  the 
right  side  of  Figure  4  corresponds  to  the  viewing  lens  of 
Figure  2  (on  the  left) . 

Selection  of  nominal  characteristics  of  various  optical 
components  used  in  the  device  is  made  during  design  of  the 
device.  Determination  of  actual  characteristics  of  these 
components  is  made  during  alignment  and  assembly  of  the  device. 
The  six  parametric  constants  determined  during  assembly  of  the 
device — aL,  aR,  ,  RR,  PL,  and  PR — are  shown  in  Figure  4. 

These  constants  are  used  to  uniquely  specify  the  locations  of 

the  various  dark  bands  produced  by  the  moire  device.  Following 
is  a  brief  description  of  each  of  these  parameters.  Subscripts 
R  and  L  refer  to  the  right  and  left  side  of  Figure  4. 


Parameter 


Description 


a 


Angle  optical  axis  makes  with  X-axis 


R  Normal  distance  from  center  of  lens  to 

X-axis  or  target  surface  (or  tangent  to 
target  surface  in  case  of  curved  targets) . 
In  actual  practice  this  value  is  calcu¬ 
lated  using  the  relationship  shown  in 
Figure  4  and  the  appropriate  values  of  a 
and  the  distance,  d,  which  was  measured 
during  alignment  of  the  individual  heads. 
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Parameter 


Description 

P  A  pseudopitch  used  when  determining,  in 

slope-intercept  form,  the  equations  of 
lines  representing  the  various  rays.  This 
dimensionless  value  is  determined,  as 
shown  in  Figure  4,  using  the  Ronchi  ruling 
pitch  (furnished  by  Ronchi  ruling  manufac¬ 
turer)  ,  the  system  magnification  (measured 
during  the  alignment  of  the  individual 
heads) ,  and  the  previously  mentioned 
distance,  d. 

Target  surfaces  shown  in  the  four  figures  just  presented 
are  represented  as  planar  surfaces  which  are  tangent  to  fringe 
zero.  This  orientation  is  not  required  for  satisfactory  opera¬ 
tion  of  the  device,  and  practically  speaking,  is  rarely  achieved. 
When  the  device  is  used  with  convex  or  concave  surfaces,  the 
preferred  orientation  of  the  device  has  the  displacement  axis 
of  Figure  3  coincide  with  a  line  normal  to  a  plane  tangent  to 
the  target  surface  at  the  point  where  intersection  of  the  X-axis 
interruptions  in  the  locator  lines  occurs.  The  relationship 
of  the  moire  device  displacement  axis  (fixed  by  optical  compo¬ 
nents  in  the  moire  device)  and  a  line  normal  to  the  target  at 
the  point  of  tangency  just  described  is  the  preferred  orienta¬ 
tion  of  the  device  and  target  for  all  applications  since  it 
facilitates  identification  of  fringe  zero.  When  the  device  is 
mounted  independent  of  the  target,  the  user  must  carefully 
position  the  device  for  each  test  to  insure  this  orientation 
is  achieved  as  closely  as  possible.  In  other  applications, 
e.g.,  measurement  of  aircraft  transparency  deflections,  the 
moire  device  may  be  installed  in  mounts  which  may  or  may  not 
provide  for  adjustment  of  the  relationship  of  the  device  to 
the  target  surface.  In  these  cases,  the  general  relationship 
of  the  target  and  the  device  is  fixed  during  design  of  the 
device  mounting  hardware . 

Most  uses  of  the  moire  device  require  some  preparation 
of  the  target  surface.  Typical  preparations  consist  of  the 
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application  of  a  reflective  coating,  usually  white  paint,  and 
the  installation  of  fiducial  marks  on  the  reflective  screen. 

The  fiducial  marks  are  used  to  provide  data  for  film  magnifica¬ 
tion  determinations.  A  suggested  fiducial  line  layout  is 
presented  in  Figure  5.  Use  of  X's  rather  than  short  line  seg¬ 
ments  as  reference  points  on  the  horizontal  axis  eliminates 
the  likelihood  of  one  or  more  of  the  line  segments  being 
obscured  by  the  lines  in  the  projected  or  viewed  Ronchi  ruling. 
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Figure  5.  Suggested  Fiducial  Pattern  for  Use  with  Moire  Device. 


SECTION  III 
DATA  REDUCTION 

Data  used  to  generate  deflection  information  are  obtained 
from  various  photographs  taken  during  the  experiments  or  tests. 

A  single  photograph  may  be  analyzed  in  the  case  of  a  statically 
loaded  target.  When  the  target  is  impulsively  loaded,  data  are 
taken  from  photographs  of  selected  frames  of  a  high-speed  film 
of  the  transient  fringe  patterns.  Measurements  of  the  locations 
of  timing  marks  placed  on  the  high-speed  films  are  used  to 
determine  the  framing  rate  of  the  camera  and  the  time  interval 
between  successive  fringe  patterns  for  these  dynamic  events. 

In  this  section,  detailed  instructions  for  making  a  variety  of 
measurements  are  provided.  The  measurements  are  described  and 
presented  in  the  order  in  which  they  are  entered  into  the  com¬ 
puter  program.  A  listing  and  description  of  sample  data  items 
is  presented  in  Section  V.  The  first  data  entries  consist  of 
specific  moire  device  parameters  and  information  unique  to  the 
test.  An  explanation  of  each  of  these  data  items  is  given 
in  Section  V. 

Although  the  original  negative  film  may  be  used  for  data 
reduction  purposes,  analysis  of  the  moire  patterns  is  greatly 
facilitated  by  using  enlarged  prints  of  the  individual  frames. 

In  the  case  of  multiple  prints,  it  is  extremely  important  that 
all  frames  taken  from  a  data  film  be  printed  at  the  identical 
magnification.  Data  used  to  obtain  film  magnification  infor¬ 
mation  are  taken  from  pre-event  photographs.  Post  event 
photographs  are  used  to  obtain  measurements  of  fringe  location 
and,  ultimately,  target  displacement.  Since  the  fiducial  lines 
are  applied  to  a  surface  which  is  in  motion  during  the  event, 
film  magnification  information  taken  from  post  event  frames 
should  be  regarded  as  suspect. 

Because  the  axis  of  the  viewing  system  is  at  an  angle  to 
the  target  surface,  distortion  of  the  image  produced  at  the 
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viewing  system  Ronchi  ruling  occurs.  Consequently,  magnifica¬ 
tion  of  the  target  image  varies  continuously  in  the  horizontal 
(X)  and  vertical  (Y)  directions.  Since  out-of-plane  displace¬ 
ment  data  are  derived  from  measurements  taken  from  the  film, 
conversion  of  measured  dimensions  (film)  to  true  dimensions 
(target)  requires  that  the  variation  in  film  magnification  be 
accounted  for  properly. 

Film  magnification  along  the  X-axis  is  determined  using 
measurements  taken  from  one  or  more  of  the  frames  immediately 
preceding  the  event.  (Only  one  set  of  averaged  values  is 
entered  into  the  program,  however.)  Profiles  other  than  those 
along  the  X-axis  are  usually  taken  from  sections  parallel  to 
the  X-axis  and  through  the  short  horizontal  lines  on  the  two 
long  vertical  lines  of  the  fiducial  system  (Figure  5) .  Mag¬ 
nification  in  the  Y  direction  is  thus  treated  graphically  and 
no  further  corrections  are  required.  An  illustration  of  a 
modified  target-viewing  head  relationship  is  presented  in 
Figure  6.  In  this  figure,  film  normally  in  the  camera  is  shown 
in  the  plane  of  the  viewing  head  Ronchi  ruling.  A  view  of  the 
undistorted  target  surface  and  fiducial  system  is  shown  in 
the  upper  right  of  this  figure.  A  view  of  the  distorted  image 
produced  on  the  film  is  shown  at  the  lower  right  of  Figure  6. 
Also  shown  in  this  figure  is  the  position  of  the  interruption 
in  the  locator  line  with  respect  to  the  X-axis. 

As  shown  in  the  lower  right  portion  of  Figure  6,  the 
locator  line  (ray  zero)  does  not  coincide  with  the  center  of 
the  fiducial  pattern  and  an  offset  of  unknown  magnitude  results. 
The  magnitude  of  this  offset  and  actual  film  magnification 
data  are  determined  graphically  using  measured  values  taken  from 
the  pre-event  photographs.  All  measurements  are  made  with 
respect  to  the  locator  line  of  the  viewing  head.  This  line 
will  always  appear  as  a  straight  line  in  all  photographs.  Raw 
film  magnification  data  consist  of  a  table  of  measured  distances 
to  known  points  (X's  on  the  horizontal  fiducial  line)  and  the 
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RAY  ZERO 


NORMAL  VIEW  OF  TARGET 
AND  FIDUCIAL  LINES 


VIEW  OF  DISTORTED  TARGET 

IMAGE  RECORDED  ON  FILM - - 

(FRINGE  PATTERN  OMITTED 
FOR  CLARITY) 


Figure  6-  illustration  oi  Technique  wnich  Produces  a  Distorted 
Image  of  Target  Surface  and  Fiducial  Pattern.  Refer¬ 
ence  to  Figures  3  and  4  will  clarify  the  apparent 
reversal  of  right  and  left  in  the  distorted  image  of 
the  target  surface  and  fiducial  pattern. 

corresponding  true  distance  to  these  points.  Measured  distances 
to  the  right  of  the  locator  line  have  positive  values  and  mea¬ 
sured  values  to  the  left  of  the  locator  line  have  negative 
values.  When  these  data  are  plotted  as  shown  in  Figure  7,  the 
magnitude  of  the  offset  between  the  center  of  the  moire  device 
and  the  center  of  the  fiducial  pattern  corresponds  to  the  true 
distance  between  the  parallel  lines  representing  the  center  of 
the  moire  device  and  the  center  of  the  fiducial  system. 

The  value  of  offset  thus  obtained  is  used  to  "correct" 
the  true  values  and  generate  a  table  of  "corrected"  film  mag¬ 
nification  data  which  are  entered  into  program  PANDIS.  For 
example,  the  true  position  of  the  second  X  to  the  right  of 
center  of  the  fiducial  pattern  is  +4  inches.  This  true  value 
is  corrected  by  subtracting  the  algebraic  value  of  the  offset 
from  the  true  position  of  the  fiducial  mark,  i.e.,  +4  -  (+1.2)  = 
+2.8.  All  other  true  values  are  corrected  in  a  like  manner. 


+  3 


Figure  7.  Illustration  of  Graphical  Method  Used  to  Determine 
Offset  and  "Corrected"  Film  Magnification  Data. 

The  measured  distance  from  the  locator  line  to  this  X  (and  all 
other  X's  in  the  photograph)  remains  unchanged. 

Processing  of  raw  film  magnification  data  in  PANDIS  yields 
the  coefficients  of  a  second  order  polynomial.  Subsequent  con¬ 
version  of  measured  distances  taken  from  data  frames  to  true 
distances  is  accomplished  by  evaluating  the  polynomial  using 
the  measured  values  as  the  independent  variable. 

The  framing  rate  of  the  high-speed  camera,  in  pictures  per 
second  (pps) ,  is  determined  for  that  portion  of  the  high-speed 
film  in  which  response  of  the  target  is  recorded.  Framing  rate 
is  determined  from  measurements  illustrated  in  Figure  8.  In 
this  figure,  timing  mark  zero  was  arbitrarily  selected  as  the 
sixth  mark  before  the  frame  in  which  response  of  the  target  was 
first  observed.  In  practice  any  of  the  several  marks  appearing 
before  this  frame  may  be  selected  as  timing  mark  zero.  All 
subsequent  measurements  are  made  with  respect  to  timing  mark 
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Figure  8.  Illustration  of  a  Portion  of  High-Speed  Film  Showing 
Measurement  Points  for  Film  Speed  Determination. 

zero.  In  program  PANDIS,  the  time  to  successive  marks,  R(  ) , 
and  the  corresponding  measured  distances,  U(  ),  are  entered 
as  data  pairs  and  used  in  conjunction  with  certain  physical 
characteristics  of  16  mm  film  to  determine  the  coefficients 
of  a  polynomial  which  yields  the  instantaneous  framing  rate 
of  the  camera  as  a  function  of  frame  number. 

The  frame  of  film  immediately  before  the  frame  in  which 
response  of  the  target  is  observed  is  identified  as  frame  zero. 
Frames  exposed  after  frame  zero  are  numbered  sequentially  in 
ascending  order.  The  first  measurement  to  be  taken,  U(0),  is 
used  to  determine  the  location  oij  frame  zero  with  respect  to 
timing  mark  zero.  The  other  measured  distances,  U(l),  U(2), 

U( 3),  etc.,  are  used  to  determine  the  actual  framing  rate  of 
the  camera.  Note  that  the  distance  to  timing  mark  zero,  U(l), 
is  zero  and  entered  accordingly  in  statement  100  of  the  Sample 
Data  (Section  V) . 


Since  the  system  of  dark  bands  generated  by  the  moire 
device  is  unique  for  the  set  of  system  parameters,  displacement 
of  a  target  at  any  point  can  be  determined  if  the  X-coordinate 
and  fringe  number  of  the  point  are  specified.  The  X-coordinate 
and  fringe  number  are  obtained  from  analysis  of  a  photograph 
of  the  fringe  pattern.  A  simulated  photograph  of  an  ideal 
target  and  fringe  pattern  is  presented  in  Figure  9.  The  intended 
use  of  the  out-of-plane  displacement  data  determines  the  number 
of  measurements  taken  from  the  photograph.  If  the  deflected 
shape  along  one  section  is  required  (usually  along  the  horizon¬ 
tal  fiducial  line) ,  measurements  are  taken  where  the  dark  bands 
intersect  the  fiducial  line.  If  a  map  of  the  deflected  surface 
is  desired,  additional  profiles  are  obtained  along  sections 
parallel  to  the  horizontal  fiducial  line.  Examples  of  both 
types  of  measurements  are  illustrated  in  Figure  9. 

Analysis  of  photographs  begins  with  identification  of 
the  various  dark  bands  in  the  moire  pattern.  Use  of  the  locator 


Figure  9.  Simulated  Photograph  of  Ideal  Moire  Pattern  Showing 
Various  Measurement  Locations. 
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lines  aids  in  the  identification  of  fringe  zero.  Study  of 
Figures  3  and  4  will  assist  in  the  identification  of  all  other 
dark  bands  in  the  pattern.  If  deflected  surface  data  are  re¬ 
quired,  lines  are  drawn  on  the  photograph  through  the  short 
horizontal  lines  which  are  part  of  the  vertical  fiducial  lines. 
Finally,  measurement  of  the  various  distances  is  made.  As  shown, 
these  distances  are  measured  from  the  center  of  the  dark  band 
to  the  locator  line.  Note  that  all  measurements  are  made  per¬ 
pendicular  to  the  locator  line  and  that  the  correct  algebraic 
sign  of  the  measured  value  must  be  observed. 

Data  for  the  measured  points  along  each  section  line  are 
entered  as  fringe  number-measured  value  pairs  beginning  with 
the  pair  at  the  left  of  the  photograph  and  ending  with  the  pair 
at  the  right  of  the  photograph.  When  data  for  a  number  of 
sections  from  a  photograph  are  to  be  entered,  data  from  the 
uppermost  section  are  entered  first;  order  of  entry  proceeds 
successively  to  the  lower  sections. 

Procedures  used  in  the  collection  of  fringe  number- 
measured  value  data  are  identical  for  the  pre-event  and  post 
event  photographs  with  one  exception.  An  additional  measurement 
is  taken  from  the  first  post  event  frame  or  the  pre-event  frame 
if  the  target  is  suitably  identified.  This  measurement,  XO  in 
Figure  9,  is  the  distance  from  the  locator  line  to  the  point 
of  the  event  or  the  center  of  the  area  of  interest.  The  value 
of  XO  is  entered  into  PANDIS  prior  to  entry  of  the  pre-event 
and  post  event  data. 

Examination  and  study  of  the  sample  data  in  Section  V 
should  clarify  and  further  the  program  user's  understanding  of 
data  entry  procedures  used  with  PANDIS. 


SECTION  IV 
PROGRAM  DESCRIPTION 


1.  PROGRAM  USAGE 

Program  PANDIS  consists  of  a  main  program  and  two  sub¬ 
routines,  LSQ  and  COR.  The  main  program  directs  reading  of 
input  variables,  raw  data,  order  of  processing  raw  data,  and 
printing  of  results.  Subroutine  LSQ  is  called  upon  several 
times  during  execution  of  the  main  program  to  fit  data  sets  to 
least-squares  polynomials.  Subroutine  COR  is  called  after  each 
raw  data  pair  has  been  read  and  processed  to  yield  a  displace¬ 
ment  value  unique  to  the  fringe  number  and  X-coordinate  of  the 
data  pair.  In  subroutine  COR,  the  initial  value  of  displacement 
is  examined  and  corrected,  if  necessary,  for  magnification 
errors  which  are  introduced  when  the  target  surface  (hence, 
measurement  plane)  has  been  displaced  a  sufficient  distance  from 
the  initial  reference  (pre-event)  plane.  A  simplified  flow 
chart  for  the  main  program  is  presented  in  Figure  10. 

Input  variables  and  test  peculiar  data  are  read  at  the 
beginning  of  the  main  program.  Certain  of  these  data  are  used 
to  characterize,  in  tabular  form,  the  location  (in  displacement, 
D,  and  X-coordinates)  of  the  fringes  which  were  observed  in  the 
test  film.  In  the  next  portion  of  the  program,  data  used  to 
determine  film  magnification  and  speed  are  processed  to  yield 
least-squares  polynomials  which  describe  the  variations,  with 
distance  and  time,  respectively,  of  these  film  parameters.  Work 
on  actual  fringe  pattern  data  begins  with  data  taken  from  the 
pre-event  frame.  Individual  sets  of  X-D  data  taken  from  this 
frame  are  fit  to  least-squares  polynomials  which  attempt  to  des¬ 
cribe  the  shape  of  the  undeflected  target  surface  along  a  partic¬ 
ular  measurement  section  line.  At  the  completion  of  printout  of 
the  results  of  the  first  fit,  the  user  is  asked  if  the  next 
higher  order  of  fit  should  be  made.  A  positive  response 
results  in  another  pass  through  subroutine  LSQ.  Printout  of 
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the  results  of  this  fit  is  again  followed  by  the  opportunity  to 
request  the  next  higher  order  of  fit.  After  various  orders  of 
fit  have  been  examined,  the  user  is  asked  to  select  that  order 
of  polynomial  which  best  describes  the  pre-event  panel  shape. 
Upon  receipt  of  a  reply,  control  of  processing  of  data  from 
post  event  frames  and  the  printing  of  results  is  returned  to 
the  main  program. 


2.  DESCRIPTION  OF  VARIABLES 

Following  is  a  listing  (and  brief  description)  of  the 
simple  variables,  arrays,  and  string  variables  used  in 
program  PANDIS.  Where  appropriate,  units  and  applicable 
reference  figures  are  included  in  this  listing. 

Simple  Variables 

A1  -  Angle,  oir,  which  optical  axis  of  R.H.  optical  system 
makes  with  X-axis,  radians  (Figure  4) . 

A2  -  Angle,  aL,  which  optical  axis  of  L.H.  optical  system 
makes  with  X-axis,  radians  (Figure  4) . 

A3  -  Angles  which  rays  in  R.H.  optical  system  make  with 
X-axis,  radians. 

A4  -  Angles  which  rays  in  L.H.  optical  system  make  with 
X-axis,  radians. 

A5  -  Angle,  a,  used  when  "correcting"  displacements  for 
magnification  errors,  radians  (Figure  11) . 

A6  -  Angle,  6,  used  when  "correcting"  displacements  for 
magnification  errors,  radians  (Figure  11). 

E  -  Standard  error  of  estimate  of  least-squares  fit. 

FO  -  Time  after  event,  evaluated  for  each  frame  number, 
msec . 

FI  -  Largest  positive  fringe  number  used  in  evaluation. 

F2  -  Largest  negative  fringe  number  used  in  evaluation. 

F3  -  Average  framing  rate  of  camera  over  an  arbitrary 

1  msec  interval,  IQ,  pps . 


18 


F4  -  Average  framing  rate  of  camera  over  the  arbitrary 
1  msec  interval,  I  +  1,  pps. 

F5  -  Aq  term  of  polynomial  describing  film  speed. 

F6  -  Ai  term  of  polynomial  describing  film  speed. 

F7  -  A2  term  of  polynomial  describing  film  speed. 

I  -  General  operator. 

II  -  Operator  to  sequence  through  all  frames  of  film  after 

pre-event  frame. 

12  -  Operator  to  sequence  through  all  sections  of  a  frame. 
J  -  General  operator. 

K  -  General  operator. 

L  -  General  operator. 

M  -  General  operator. 

Ml  -  Order  of  fit  desired  in  least-squares  subroutine. 

M2  -  Ml  +  1. 

M3  -  Aq  term  of  polynomial  describing  film  magnification, 
inches . 

M4  -  A^  term  of  polynomial  describing  film  magnification. 

M5  -  A2  term  of  polynomial  describing  film  magnification, 

inches”! . 

N  -  General  operator. 

NO  -  Number  of  frames  of  film  to  be  evaluated  (including 
pre-event  frame. 

N1  -  N9  +  1. 

N2  -  Frame  number. 

N3  -  Number  of  sections  in  a  frame. 

N4  -  Value  used  to  signal  that  insufficient  data  are 
available  for  a  least-squares  fit. 

N9  -  Number  of  rays  on  either  side  of  optical  axis 
(Figure  4) . 
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PO  -  Proportionality  factor  used  during  interpolation  in 
X-D  Tables. 

PI  -  Number  of  data  pairs  for  a  particular  use  of  least- 
squares  subroutine  or  other  functions  requiring 
knowledge  of  number  of  data  pairs  for  a  particular 
operation . 

P7  -  Page  counter. 

P9  -  Pi  (3.14159) . 

R1  -  Normal  distance  from  R.H.  optical  system  lens  to 
X-axis,  mm  (Figure  4). 

R2  -  Normal  distance  from  L.H.  optical  system  lens  to 
X-axis,  mm  (Figure  4). 

R3  -  Constant  used  when  "correcting"  displacements  for 
magnification  errors,  mm. 

R8  -  Pseudopitch  of  R.H.  optical  system  (Figure  4). 

R9  -  Pseudopitch  of  L.H.  optical  system  (Figure  4). 

S  -  Summation  in  least-squares  subroutine. 

SO  -  Calculated  displacement  (w/r  to  Moire  Device)  using 
pre-impact  data  fit,  mm. 

S2  -  Ax  in  "correction"  of  displacements  for  magnification 
errors,  mm. 

56  -  Temporary  storage  of  X  values,  mm. 

57  -  Temporary  storage  of  D  values,  mm. 

XO  -  Point  of  impact  w/r  to  centerline  of  Moire  Device, 
inches  and  mm  (Figure  9) . 

YO  -  Counter  in  correction  subroutine. 

Y2  -  Summation  in  least-squares  subroutine. 

Y8  -  Projectile  weight,  variable. 

Y9  -  Impact  velocity,  variable. 

Z  -  Absolute  value  of  (Y2-S) . 
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F  (  ) 

G  (  ) 


P  (  ) 


Q  (  ) 


R  (  ) 


-  Counter  number  for  X  and  D  values  tabulated 
for  each  fringe. 


-  Storage  for  various  summations  in  least- 
squares  subroutine. 


H  (  ,  ) 

M  (  ) 

N  (  ) 

0  (  ,  ) 


-  Specific  fringe  numbers  from  pre-event  frames. 

-  Slopes  of  rays  in  right  optical  system. 

-  Slopes  of  rays  in  left  optical  system. 


Storage  of  coefficients  of  polynomials  des¬ 
cribing  each  section  in  pre-event  frame. 
0(9,  )  =  PI  and  0(10,  )  =  Ml. 


-  X-values  in  least-squares  subroutine,  inches 
or  mm;  measured  values  for  film  magnification 
data,  inches;  and  measured  values  for  post 
event  frames,  inches. 


Y-values  in  least-squares  subroutine,  inches 
or  mm;  true  values  for  film  magnification 
data,  inches;  and  fringe  numbers  for  post 
event  frames. 


Time  values  for  film  speed  determination, 
msec  and  measured  value  for  pre-event  data, 
inches . 


T  (  ,  ) 

U  (  ) 


-  Various  Exi^Vj.  in  least-squares  subroutine. 


Measured  distances  for  film  speed  determination 
inches  (Figure  8) ,  and  fringe  number  for  pre¬ 
event  data. 


V  (  ,  )  -  Inverse  of  Y( 


)  in  least-squares  subroutine 


W  (  ,  )  -  Specific  X-values  from  pre-event  frame,  mm. 

X  (  ,  )  -  X-values  in  table  of  fringe  coordinates,  mm. 

Y  (  ,  )  -  Various  £x^m  in  least-squares  subroutine. 

String  Variables 

D$  -  Date. 

S$  -  Shot  or  test  identification  code. 

X$  -  Query  for  next  order  of  fit  in  least-squares 

fit  for  pre-event  data. 

P$  (  )  -  P$(0)  -  Target  description. 

P$(l)  to  P$ (9)  -  ST,  ND,  RD,  TH,  TH,  TH,  TH, 
TH,  TH  suffixes  added  to  order  of  fit  to  yield 
1st,  2nd,  etc. 

P$(10)  -  Projectile  description. 

P$(ll)  to  P$(20)  -  Section  location 
designators . 

3.  PROGRAM  LISTING 

Following  is  a  listing  of  the  BASIC  statements  which 
comprise  program  PANDIS .  The  program,  as  listed,  was  run  on 
a  UNIVAC  90/80-3  computer  using  the  sample  data  included  in 
the  next  section  to  produce  the  printout  presented  at  the  end 
of  the  next  section.  The  listing  has  been  separated  into  a 
number  of  statement  groups.  In  general,  each  group  corresponds 
to  one  of  the  operations  shown  on  the  flowchart  in  Figure  10. 
Where  appropriate,  a  short  explanation  of  the  function  of  each 
group  of  statements  is  presented  prior  to  the  listing  of 
statements . 


0  REM  STATEMENT  NUMBERS  1  THROUGH  999  ARE  AVAILABLE  LOR  DATA 


See  Sample  Data  (Section  V)  for  format  requirements. 
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The  following  statements  allocate  storage  for  the  variously 
named  arrays . 


1000  tun  Ado,! ) »  p  < 200) ,  c  < :.’oo >  ,d<  1.00,30)  ,f<3<>,  io> 

(010  II  r H  F  <200) ,G< 10) ,H( 30, 10) ,H<200> ,N< 200) ,0(10,10) 

1  020  IUM  50)  .|-1\  20  )  ,0(50)  ,R<50>  ,  T>  10,1  ) 

10.50  tun  11(50)  .0(10,10)  ,U(30. 10)  ,  >.  <  100,30)  ,  Y<  10, 10) 

2000  RFM  PARAMETRIC*  CONSTANT:;  FOR  MOIRE  APPAKAlT  10 

2  010  R  FI  A  It  R 1  ,  R  2 , 0 9 ,  R  8 ,  R  V 
2020  READ  A1 ,A2,F1 • F2 . N3 

20.50  PV-,5, 1  4159 
2040  N1  - N°  t  I 
2050  F  OR  I  -1  TO  9 
20 AO  RE  Alt  F  t  <  F  ) 

2070  NF  X  I  I 

2000  R3=R1 /TAN  <  A  I > 

2090  FOR  1  =  1  I  10  M3 1  10 
2100  RE  AD  F't  (  1 ) 

2110  NEXT  1 

2200  RFM  -  READ  ill  .Ml  RAI  INFORMATION  REGARDING  SHOT 
2210  RE  AP  S*»D*»P*<0> 

2220  RF Alt  P$< 10) »YH»Y9»N0 

In  the  next  operation,  the  general  form  of  equation  used 
to  determine  the  slopes  of  the  various  rays  is: 


m  =  tan  6 

where  m  is  the  slope  of  the  ray  and  6  is  the  angle  which  the  ray 
makes  with  the  X  axis.  The  D-intercept,  b,  of  each  ray  is 
computed  using  the  following  relationship  and  values  illustrated 
in  Figure  4. 


In  the  program,  care  is  taken  to  insure  that  the  proper  algebraic 
signs  are  observed  during  computation  of  the  various  values. 


2400  REM  -  DETERMINE:  SLOPES  AND  INTERCEPTS  OF  ALL  RAYS 

2410  FOR  1=  N9  TO  NV 

2420  .1=1  TNI 

7430  A3=ATN< I*R8> 

2440  A4=ATN<T*R9> 

2450  M< J) -TAN ( At  A3 ) 

24A0  »<J)=R1  <  (R.1 /IANFA1  )  )*M<  J)  ) 

2470  N ( J )  =--TAN<  A2 1 A4  ) 

2480  C  (  J )  =R2+  <  < R2/TAN  <  A2  >  >  *N<  .1 )  ) 

2490  NEXT  I 


Before  the  coordinates  of  the  intersections  of  rays  are 
tabulated,  the  fringe  number  of  the  intersection  is  checked  to 
determine  whether  or  not  it  falls  within  the  range  of  fringe 
numbers  that  appeared  on  the  data  photographs.  Intersections  for 
fringe  numbers  which  lie  outside  the  range  are  not  computed. 
Computation  of  the  fringe  numbers  and  X  and  D  coordinates  of 
the  various  intersections  is  made  using  relationships  given  in 
Figure  4. 
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polynomial  in  subroutine  LSQ.  The  coefficients  of  the  resulting 
polynomial  are  printed  and  stored  for  later  use. 
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4000  REM  -  REAP  Ell  M  MAGNIFICATION  DATA  AND  F  IT  TO  A  2ND-  ORDER  POLYNOMIAL 
4010  REAP  PI, Ml 
4020  FOR  1=1  TO  PI 
4030  REAP  r<T>'0(I> 

4040  NEXT  I 
4050  COSUP  10000 
4060  M3  =  A  < 1 , 1 > 

4070  M4=A (2,1) 

4080  M5=A<  3 , 1 ) 

4090  LIST  USING  91 70, M3, M4 ,M5,E 


In  statements  4450  and  4480,  measured  values  taken  from  the 
film  (Figure  8)  are  used  to  compute  "instantaneous”  framing 
rates,  F3  and  F4 ,  for  successive  time  intervals.  These  computa¬ 
tions  are  tailored  to  16  mm  film  (40  frames  per  foot)  with 
timing  marks  applied  to  the  film  at  1  msec  intervals.  Successive 
"instantaneous"  framing  rates  are  averaged  to  provide  the  framing 
rate  data  which  are  fit  to  a  polynomial  in  subroutine  LSQ. 
Coefficients  of  the  polynomial  are  printed  and  stored  for  later 
use. 


4400  REM  REAP  FILM  SPEED  DATA  AND  FIT  TO  A  2ND  ORDER  POLYNOMIAL 
4410  READ  P1»M1»U<0> 

4420  FOR  1=1  TO  PI 
4430  READ  R(I),U<I) 

4440  NEXT  I 

4450  F3=<  <U<2)-U< 1 ) ) *10000 ) / < 3* ( R < 2 )  R ( 1 ) ) > 

4460  FOR  1=1  TO  PI -2 

4470  F*  <  I  >  ==  <  LI  <  I  *  1 )  IK  0 )  )  #10/3 

4400  F4=<  <U< 1+2)  U ( I  + 1 ) >*1 000 0 ) / < 3* (R ( I F2 ) -R < I  +  1 ) )  ) 

4490  0(I)=(F3+F4)/2 
4500  F3=F4 
4510  NEXT  I 
4520  Plapl-2 
4530  GOSUB  10000 
4540  F5=  ACI.  ,1) 

4550  F6=A (2,1 ) 

4560  F7=A<3» 1 > 

4570  LIST  USING  9 180 , F5, F6 , F7 , E 


4800  REM  REAP,  COMPUTE.  AND  PRINT  REAL  VALUE  OF  POINT  OF  IMPACT 
4810  READ  XO 

4820  X0= ( M3 1 M4*X0+M5*X0**2 ) *25 . 4 


4830  LIST  USING  9670 
4840  IF  SON ( XO )  0  THEN  4870 
4850  LIST  USING  9190«ADS(X0> 
4860  GOTO  4880 


4870  LIST 
4800  LIST 
4890  LIST 
4900  LIST 
4910  LIST 


USING  9700. APS fXO) 
USING  9260 
USING  9270 
USING  9280 
USING  9290 


4920  I.  1ST  USING  9680 
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Raw  data  for  the  pre-event  frame  are  processed  on  a 
section-by-section,  point-by-point  basis.  After  each  fringe 
number-measured  value  pair  has  been  read,  the  measured  value 
is  processed  to  compensate  for  film  magnification  errors  and 
to  convert  units  of  measurement.  Using  the  fringe  number  and 
X-coordinate  as  access  controls,  a  displacement  value  is  ob¬ 
tained  from  a  linear  interpolation  of  tabulated  fringe  coor¬ 
dinates.  These  initial  values  of  X  and  D  are  stored  and  then 
processed  in  subroutine  COR  to  determine  the  magnitude  of 
correction,  if  any,  which  must  be  applied  to  the  raw  values. 
After  the  corrections  for  film  magnification  have  been  made, 
the  raw  and  corrected  values  (and  a  count  of  the  number  of 
adjustments  made  in  subroutine  COR)  are  printed  at  the  terminal. 


'■.ooo  Kin  id  on  Hirwi.r  on.  ,m  ,,-.md  n  ■  0010  nu<  i  id:  iriivu  i  rord  i  n..  10 
sol.,  i  iik  r,.1  i  in  d > 

*.u,'0  I 'kid  i 

'.o  r-k i  in  l  a  (  i/i  i  o  i 
soio  i  k in i 

SOSO  I  I  I  rJ  I  1 1'.  |  ill.  74  10 
SO..')  I'kJ  NT  If.,  I  Nil  v  I  ’  , 

SO'O  PRINT 
'.OHO  k'l  oil  I  I 
S0"0  I  III-  I  I  III  I'  I 
5100  Id  On  mi.  .I.'ii  '. 

SI  io  k<  I  )- «  n  vih -Ukf  l  >  ! h5*k  I  )  \>>*r>5.4 
s  I  .'<)  so  -!•<  i  > 

S  I  ,)0  .1=1  I  11(1)11 

SI  40  I  Ok  k  I  III  N'MI 

SI  so  If  X<l... I)  F'<  I  )  I  III  N  5170 

5  I  AO  NEXT  k 

s  I  /'()  PO  (  X  !  I.  I  r  J )  r  (  I  >  >  /  I  X  Ik  1.1)  X  (  l\  r  J  )  I 
'll  UO  0CI>  HO,  I  -  .1)  1 1  <  II  (k  ,  I .  Otk  I  •  I  >  > 

s  i'.'o  D 

5/00  ODSIIIi  I  5000 
5710  k  <  I  )  :=  p  (  T  ) 

si>  ’o  u<  i ,  i )  in  i ) 
s;mo  f<  i  i  on  i; . 

5.’ 40  in  J  .  17  .  I  (  I  ) 

■..'so  ii(  i;  .=()■  r  i 
5 7 AO  fit. r  I  )~-U( i ) 

57/0  PRINT  USING  9700 1 ) 1 ) »S7»U<  I > .Yd 
5 7 HO  Nf  X  I  I 


The  corrected  X  and  D  data  pairs  for  each  pre-event 
section  are  fit  to  a  second  order  polynomial  in  subroutine 
LSQ .  Coefficients  of  the  resulting  polynomial  as  well  as  the 
measured  and  fitted  data  are  printed  at  the  terminal  and  the  user 


is  asked  if  the  next  higher  order  of  fit  is  desired.  A  positive 
response  results  in  a  repetition  (at  the  next  higher  order)  of 
the  curve  fitting  and  printing  process.  A  negative  response 
terminates  the  procedure  and  the  user  is  asked  to  specify  the 
order  of  fit  which  will  be  used  to  describe  the  particular 
section.  Pre-event  displacement  data  are  fit  to  a  polynomial  to 
provide  a  continuous  functional  relationship  which  is  used  during 
determinations  of  the  relative  displacement  of  the  panel  for  each 
of  the  varying  X  values  from  the  post  event  data. 


5400  REM  -  FIT  X.D  DATA  TO  LEAST  SQUARES  POLYNOMIALS  AND  PRINT  RESULTS 

5410  Ml =2 

5420  N4=0 

5430  G0SUR  10000 

5440  IF  N4=-1  THEN  5710 

5450  PRINT 

5460  PRINT 

5470  PRINT  USING  9410.MI  .PKM1  ) 

5480  PRINT  USING  9220 

5490  FOR  J- 0  TO  Ml 

5500  PR  I  NT  US  I NG  92.10  .  J  t  A  (  J  *  t.l> 

5510  NEXT  J 

5520  PRINT  USING  9240. E 
5510  PRINT 

5540  PRINT  USING  9440 
5550  PRINT  USING  9450 
5560  PRINT 
5570  FOR  J=1  TO  PI 
5580  SO  "A  <1  » 1  > 

5590  FOR  K  -I  TO  Ml 

5600  SO-- SO  I A  ( K  H  .  1  )  *  (  R  <  DM 

5610  NEXT  K 

5620  PRINT  US  ING  °4AO . F  ( . |  > » R  i  I )  .  II ( .1  >  •  SO » ll<  J )  SO 
5610  NEXT  I 
5640  PRINT 
5650  PRINT 

5660  PRINT  TIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL*  (YCS/NO)*. 

56/0  INPUT  X4> 

5680  IF  X V-  " NO ■  THEN  5710 
5690  Ml -Ml 11 
5700  GOTO  5410 
5710  PRINT 

5720  PRINT  •DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE  * 

5710  PRINT  ’PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  'i 

5740  INPUT  Ml 

5750  PRINT 

5760  PRINT 

57  70  GOSIJD  1.0000 

5780  REM  STORE  COEFFICIENTS  OF  CURVE  FIT  FOR  EACH  SECTION  LINE 
5790  FOR  1=0  TO  Ml 
5800  0<  HI .  I2)=A(  1  +  1 . 1  > 

5810  NEXI  I 

5820  LIS  I  USING  9100. 12.0'  1.121.0(2. 12) . 0 ( 1 . 12 ) . 0 ( 4 . 12 ) .0(5. 12)  .E 
5810  0(10.12) =M1 
5840  0(9»I2)=P1 
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58  so  ni  a  i  i  :> 

5SV.0  RIM  COrlPI  I  ft  PR  INI  Oil  I  Of  Mill  I'M. I 

5870  ITS  I  US  INI,  '•<  '50. FI  Kill..  H.il  <  U  • 


6000 

RFM 

PRINTOUT  FIR 

ST 

PAGE  OF  DATA  <  E  OR  PUT 

IMPACT  frame: 

6010 

E  ‘  7  - 
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1. 1 S 

I 

USING 
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'.NOLI 
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1.  IS 

i 
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In  the  next  segment  of  the  program,  fringe  number- 
measured  value  pairs  for  all  post  event  frames  are  processed. 
Page  number,  frame  number,  shot  number,  and  the  time  after  the 
event  and  other  page  heading  information  are  printed  first. 
Processed  data  for  successive  sections  of  each  frame  are  then 
listed.  As  for  the  pre-event  case,  raw  measured  values  are 
corrected  for  film  magnification  before  they  are  listed.  Only 
the  corrected  values  are  listed,  however.  In  addition  to 
fringe  number  and  X  and  D  values  (in  device  coordinates)  the 
listing  includes  X  values  referenced  to  the  point  of  the 
event  and  displacement  values  referenced  to  the  pre-event 
position  of  the  target. 


7000  RE  M  RIAD  ANH  PRINT  All  FRAMES  ill  U,il  A  hI  ICR  IRE  IMPACT  FRAME 
7010  FOR  II- I  TO  NO  I 
7020  RE  AP  n:> 

70,10  pr  r  n  i  n.>; 

7040  I  1ST  USING  9  A 1 0  .  S  i  »  E '  7  »  NO  I  I 
7050  e  o  •  n:>/ < rsi r  i F7*N2"? iui ooo 

70 A0  LIST  USING  9620 • N7 - 1  0 
7070  LIST  USING  9630 
7080  LIST  USING  9640 
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7090  LIST  USING  9650 
7100  LIST  USING  Vf-.'.O 
7  J  10  FOR  1 2- I  i  n  ll  < 

7 120  Rf AM  FI 

7130  If  fl-O  WIN  73  So 

7140  LIST  USING  9710»P* < 17+ 10) 

7150  TOR  1=1  TO  FI 
71  AO  Rf  Oil  f) <  I  >  *  F  <  1  ) 

7  1  7< •  F  <  T  >  0(  |  t 

71  HO  F <  I  )  ifUlriHLi  I  '  »N5*F<  I  >"  2)*25.  4 

7  170  J-F  1  0 <  DM 

7 '700  FOR  K  -t  TO  07 1  I 

7?lo  IF  X(K.J)  F(T)  THEN  7230 

722 0  NEXT  K 

7730  PO=  <  X  <  K  If  J)  P  <  I  )  >  /  <  X  <  K  1,1/  X  <  I-  •  J  >  > 

72 AO  i:+  *  I  )  - n < K  1  »  I  <  IPG#  (  n ( K  »  J  >  IKK  I  -  I >  ' 

7750  GOSUR  15000 
7 7 AO  50=0 < 1 • [7> 

7770  FOR  I  =1  TO  o< JO, f?> 

7280  SO=SOfOfl.  1  1 «  1  ,’>»<('(]  )**1  < 

729 0  NEXT  l 

7300  LIST  USING  9A?0»I  1  I  ) «F<  1  ) »Q< I > »P* I >-XOrQ< I >~S0 
7310  NEXT  I 
7370  GOTO  7340 

7330  LIST  USING  9720 »P*< 17110) 

7340  NEXT  17 
7350  F7-F7U 

7360  RE M  GO  10  Niff  f  RArlf  hND  Fh Of  Of  PRINTOUT 

7370  Nf  X  I  I  I 

8000  FRINI  *  JOB  COMF1  f  if  El  GOING  TO  t  N»  STATEMENT 

00 10  GOTO  20000 


Following  is  a  listing  of  the  various  image  statements 
used  in  conjunction  with  the  various  PRINT  USING  and  LIST  USING 
commands.  Content  of  these  statements  can  be  changed  to  suit 
the  heed  of  the  user.  Image  statements  listed  and  used  for  the 
sample  problem  were  used  for  an  aircraft  transparency  deflection 
measurement  application. 


9000 

9010 

9020 

9030 

9040 

9050 

9060 

9070 

9080 

9090 

9100 

9110 


REM 

SI 

SO 

SO 


FOl. 

OF 


IMAGE  STATEMENTS  FOR  TITLE  PAGE  PRINTOUT 

F16-A  BIRDSTRIKE  PROGRAM 

♦#***♦♦*♦*****♦♦**#*#*»*******♦***# 
*•**»*•*•*#••*•••*•**•*•••*•• 
******************************************** 
************************************** 

*.**  LBS 
***  KNOTS 

OF  FILM  < AFTER  IMPACT)  EVALUATED  FOR  THIS  SHOT 


SHOT  NUMBER 
DATE  FIRED 
TRANSPARENCY 
PROJECTILE 
WEIGHT 
VELOCITY 
**  FRAMES 


LOWING  ARE  THE  COEFFICIENTS  OF  LEAST- SQUARES  POLYNOMIALS 
THE  FORMS 

Y  =  AO  +  A1*X  +  ...  +  AN*X#*N» 
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■  9120 


FIT  TO  THE  EXPERIMENTAL 


9130 
9  MO 
9150 
9160 
91 .0 
91(10 
9190 
9200 
9210 
9220 
9230 
9240 
9250 
9260 
9270 


5  DATA.  FILM  MAGNIFICATION ,  IN  REAL  INCHES ,  IS  GIVEN  AS  A  FUNCTION 
5  OF  MEASURED  INCHES  AND  FILM  SPEED*  IN  PICTURES  PER  SECOND*  IS 
J  GIVEN  AS  A  FUNCTION  OF  FRAME  NUMBER. 

{©COEFFICIENTS:  AO  A1  A2  STD  ER  OF  ESI 

5  0F.TLM  MAGN1FJ  *.*#*#•'!  M  --#.#*###'('  1  #.**###<!  !'  #.#**#*'!!' 

iOFIlM  SPEED!  • . ***** ■ ■ t •  -#.#####!'!'  #.###*#!!!(  #.##*##!!'' 

:  FUS  ST  A  140  LOCATED  ###.##  MM  TO  THE  RIGHT  OF  CENTER  OF  MOIRE  DEVICE 
I  I  US  SI  A  MO  LOCATED  ###,##  MM  TO  THE  LEFT  OF  CENTER  OF  MOIRE  DEVICE 
:  TNSUFFTCIFNl  DATA  FOR  FIT  TO  A  *  ##  ORDER  POLYNOMIAL 
5  PANE!  DATA  1  REFERENCED  TO  THE  MOIRE  DEVICE). 

I  Alt;  *  I  »#*»#' i'1 

I  STD.  ERROR  OF  ESI  (MATE  *. Itttt  '  <  '  ' 

{(>  PROGRAM  RUN  AT  **»**#  HRS  ON  #**##*. 

JO  **  RESULTS  OF  CURVE  FITS  TO  PRE-IMPACT  P.iNI  I  SECTIONS  *# 

;  < REFERENCED  TO  THE  MOIRE  DEVICE) 


9200  ! OSEC 
9290  ;  NO 


AID) 


COEFFICIENTS  Of 
All.) 


LEAST- SQUARE  POI  YNOHIALS 
A >.  2  >  A  i  )  A I  4  > 


STD  ERROR 
OF  ESI 


9300 
9400 
9410 
9420 
9431) 
94  40 
9  450 
9460 
9470 
9600 
96  I  0 
9690 
96.  U) 
9640 
9650 
9660 
96  70 
9600 
9690 


:  *  *.#**11(1  *.***1111  #.*##1111  #.#**1111  #.#*#'11!  -•*.*#!  !  i  I 

REM  IMAGE  STATEMENTS  FOR  PR  IN  TOUT  01  PRE  IMPACT  PANEL.  DATA 
{  FRINGE  MEASURED  X,MM  D l SI  LACK  MENT  »  MM  NO  TIMES  IN 

:  NUMBER  RAW  (  ORR  RAW  COM-  SUBR  COR 

!  COP FF IC TEN  IS  OF  A  *  #1  ORDER  POLYNOMIAL  FIT  TO  HIE  PRE- IMPACT 
!  FRINGE  MEASURED  ACTUAL  CAL CULATLD  DIFFERENCE 

I  NUMBER  X  *  MM  DIGIT  .MM  D I  SIT  *  MM  ATI  CAIX 

I  #*  ##.#♦#! 1 ! !  ##.###*' ' M  ##.*#**IMI  **.****>Mi 

! OF RAMI  NUMBER  PRE  IMPACT 

REM  IMAGE  ST  ATEMENTS  I  OR  ROUT  INE  DATA  PAGES 

: I  SHOT  NUMBER  I####*###########*##**#####**#***#  PAGE  ##  OF  #*  PAGE 


1 0  E  R  AMP  NUMBE'R  *  *  >  *******  M  S I .  C  A I  T  P I  > 

I  X  COORDINATE,  DTSPLACI  Ml  NT » 


IMPACT) 

X  COORD T WATT  - 


DISPLACEMENT* 


PRINGt  IN  MM. 

NUMBER  W/R  TO  MOIRE 
DEVICE 


IN  MM, 

W/R  TO  MOIRt 
DEVICE 


I  IT  Mh , 

W/R  TO 
EOS  ST A  140 


TN  MM, 

W/R  TO  PRE SHOT 
POSIT  ION 


#*  *#*.**« 


##.##* 


**#.#** 


**.#** 


9/00 

9/10 

9/20 


##  *11.*#  -#*#.**  *##.**  *##.**  *# 
o###*#t**#*****##*#####*##*#*##*###*##***#*#***##*##*#*##***t#*##t 
ODATA  NOT  AVAILABLE  FOR  ####*#***######*******#*#**####*###**#**###** 


In  the  following  subroutine,  LSQ,  pairs  of  data  are  fit 
to  a  polynomial  using  a  least-squares  curve  fitting  routine. 

The  sets  of  X.^  and  data  are  used  to  form  the  following  series 
of  equations: 
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EYi 

=  An 
o 

+  A. EX. 

1  l 

+  .  .  .  +  A  EX? 

m  i 

EXiYi 

=  A  EX. 

O  1 

+  Vxi 

+  .  .  .  +  A  EX?+1 

m  i 

ex1?*. 

i  i 


A  EX? 
o  1 


A. EXm+1 

1  l 


+  AmEX/ 

m  i 


In  these  equations,  n  is  the  number  of  data  pairs,  m  is  the 
order  of  the  desired  fit,  and  A0,  At,  A2  ...  are  coefficients 
of  various  terms  in  the  equations.  Since  these  equations  are 
of  the  general  form, 

y  =  A0  +  AjX  +  A2X2  +  .  .  . 

they  may  be  solved  simultaneously  to  yield  the  various  coeffi¬ 
cients.  In  subroutine  LSQ,  this  is  accomplished  by  forming 
the  various  summations,  placing  the  appropriate  values  in  arrays, 
and  using  matrix  algebra  to  solve  for  the  coefficients.  The 
various  EX?  terms  are  formed  in  statement  10140  and  later 
stored  in  array  Y;  the  various  EX^V  terms  are  formed  in  state¬ 
ment  10180  and  stored  in  array  T.  Solution  of  the  coefficients 
is  accomplished  in  statements  10350  and  10360.  Upon  completion 
of  a  determination  of  the  standard  error  of  estimate  for  the 
fit  using  the  following  relationship, 

EY2  -  EA  x'V 

Std  Error  of  Estimate  =  - — - 

n 

control  is  returned  to  the  main  program. 


ioooo  rf.m  -  least  squares  subroutine 

10010  REM  -  UP  TO  50  DATA  POINTS  MAY  BE  PIT  TO  A  9TH  ORDER  POLYNOMIAL 

10020  M2=M1 +1 

10030  IF  M2>P1  THEN  10430 

10040  MAT  A='ZER  <  M2  *  1 ) 

10050  MAT  T=ZER<M2»1> 

10060  MAT  V=ZER  <  M2  »M2 ) 

10070  MAT  Y=ZER <  M2  *  M2 ) 

10080  FOR  J=0  TO  2*Mt 


10090 

G <  J ) =0 

10100 

NEXT  J 

10110 

S=Y2=0 

10120 

FOR  N” 1  TO 

Ft 

10150 

FOR  K  -  1  TO 

2*  Ml 

10140 

G<K>  -G(M  TP 

(NV'K 

10150 

NEXT  f\ 

10140 

FOR  L  -1  TO 

M2 

10170 

IF  E  1  THEN 

10700 

10190 

ni  *  i  >=•  T<i.» 

.1  t  IliNI  *  ( 1  •  iNl'il 

10190 

00! 0  10210 

10200 

TO  •  1  >  -  1(1  r 

1  1  ICKM) 

102(0 

NEXT  1 

1 0  790 

T7  Y7  t  U  ( IJ ) ' 

*  > 

1  027,0 

n  ( o )  -  n 

10240 

NEXT  N 

1  0250 

1  M  1 

1  0  740 

r  ok  i-t  to 

M2 

10790 

i  =i.  1 1 

10280 

j-i-  i 

10290 

TOR  (.  1  TO 

1 

1  0700 

M=K-  1 

1 0  5  1  0 

N  1  .1 

1  0  >7" 

;'t  1  .  N  >  M> 

1  0,5.50 

HI  X  1  K 

1  0  540 

NI  V  T  1 

10.7,0 

MAT  '/  I NV ( r 

) 

10  >40 

MAI  A  0*1 

10.570 

TOR  [  :l  HI 

M2 

10,580 

£»•-<;»  <  A  (  T  *  1  1 

* !  :  l .  I  I  ■ 

10  5  90 

NEXT  1 

10  400 

7 '"AITS  (  Y2  0  ) 

10410 

E=S0R» 7/0(0 

• ) 

10420 

(ill  TO  10450 

104:50 

PRINT  USING 

92 10. Ml  •  I't  <  Ml  > 

10  440 

N  4  1 

10450 

RE  TURN 

Film  magnification  correction  terms  were  determined  using 
data  taken  from  the  pre-event  frame.  When  the  target  is  de¬ 
flected,  a  change  in  the  magnification  of  its  image  occurs. 

Since  the  dark  band  used  to  define  a  measurement  point  is  formed 
on  the  surface  of  the  target,  evaluation  of  the  film  magnifica¬ 
tion  polynomial  using  post  event  measurements  may  yield  results 
which  are  somewhat  in  error.  Displacement  values  obtained  with 
use  of  these  measured  values  will  also  be  in  error.  The  mag¬ 
nitude  of  the  error  increases  as  the  post  event  measurement 
plane  is  displaced  in  either  direction  from  its  pre-event 
position. 

The  error  just  described  is  reduced  to  a  specified  mini¬ 
mum  value  in  subroutine  COR  using  an  iterative  correction 
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procedure.  A  graphical  illustration  of  one  step  of  the  cor¬ 
rection  procedure  is  presented  in  Figure  11.  In  this  procedure, 
use  is  made  of  the  fact  that  each  measured  value  lies  on  a  line 
which  makes  a  unique  angle  to  the  initial  position  of  the  target 
(or  the  X-axis).  With  respect  to  the  illustration  in  Figure  11, 
angle  a  should  equal  angle  6  for  this  condition  to  be  fulfilled. 
Angle  a  (A5,  statement  15020)  is  the  true  angle  which  the  line 
of  sight  makes  with  the  pre-event  target  surface.  Angle  6  (A6, 


Figure  11.  Illustration  of  Procedure  Used  to  "Correct"  Measured 
Values  for  Magnification  Errors. 


statement  15030)  is  the  angle  which  an  apparent  line  of  sight 
makes  with  the  pre-event  target  surface.  These  angles  are 
determined  in  subroutine  COR  using  device  constants  provided 
at  the  beginning  of  the  program  and  the  particular  XA  and  D A 
values  obtained  from  the  initial  interpolation  of  the  tabulated 
X  and  D  coordinates.  The  difference  between  these  angles  is 
compared  with  a  specified  minimum  acceptable  error  (statement 
15040) .  If  the  difference  is  less  than  the  acceptable  error, 
control  is  returned  to  the  main  program. 

When  a  correction  must  be  applied,  the  magnitude  of  the 
correction,  AX  ( S 2 )  ,  is  determined  by  making  use  of  geometric 
relationships  illustrated  at  the  upper  left  of  Figure  11  and 
expressed  in  statement  15050.  The  correction  is  applied  to  the 
old  value  of  XA  and  the  new  XA  is  used  to  determine  a  new  value 
of  Da  by  interpolation  in  the  table  of  fringe  coordinates. 

Note  that  use  is  made  of  the  fact  that  the  measured  point  must 
lie  in  the  center  of  the  dark  baiid  specified  as  the  fringe 
number  in  the  original  data  pair.  The  new  values  of  X  and  D 
are  used  to  compute  a  new  angle  6.  The  difference  between  a 
and  the  new  3  is  compared  with  the  minimum  acceptable  difference 
and,  if  necessary,  the  correction  process  is  repeated.  The 
examination  and  correction  procedure  is  repeated  until  the 
difference  between  these  angles  is  less  than  the  minimum 
acceptable  value. 

15000  REM  SUBROUTINE  FOR  CORRECTION  OF  ERROR  HUE  TO  MAONIT CATION 
15010  Y0=0 

1 5020  A5---ATN  <  R 1  /  <  R3--R <  I  >  >  ) 

15030  A,  ATN<  <R1  Q<  I )  )/(R3~f'(  I )  )  ) 

15040  IF  ADS ( A5-  AA ) <5E -04  THEN  15140 

15050  S2=< < ( R3-P < I ) ) /COS <  AA ) )*TAN(A5  AA) )*SIN<A6> 

150 AO  F'<  I  )--F<  I  >+S2 

15070  FOR  K~ 1  TO  N911 

15000  IF  X!K»,JXP(I)  THEN  15100 

15090  NEXT  K 

15100  F'O-  (XCK  1 ,  J)~PO  )  )/<X(K  lt..l)~X(  KrJ)) 

15110  Q<  I  1  ,  JHFO*(D(Kf  ,.l>  -TKK  1  f  J)  ) 

15120  Y0-Y0+1 
15130  GOTO  15030 
15140  RETURN 


20000  END 


SECTION  V 

TEST  PROGRAM  AND  SAMPLE  PRINTOUT 


Except  for  order-of-fit  data  which  are  entered  from  the 
time-sharing  terminal  during  a  run,  all  data  used  during  the 
execution  of  PANDIS  are  read  from  DATA  statements.  Statement 
numbers  1  through  999  are  available  for  use  as  DATA  statements 
Following  is  a  listing  of  sample  data  which  are  read  from  these 
DATA  statements.  The  data  shown  in  these  statements  are  used 
to  produce  the  sample  printout  presented  in  this  section. 


An  explanation  of  the  data,  by  variable  name,  follows  this 
listing . 

to  I* a ( a  i v<>  jV'>.  i.;;. v.  c.".. u  ..  r-.Jtt  <■? 

.**>  DAIA  I  '  >.  i  o  ',  . 

Vi  DA  i  A  III.  III.  ill-  (!l.  Ill-  HI 

}J  DA  I A  sti  I  mil  I  I|J|  ||  AIII'.'I  liil'.v;r  1-nirJI 

32  DATA  SU  I  [fill  THI.'IHII.M  lill'Ai  l  l  iilHf 

33  DATA  r.nriION  I  IfJill  HI  I  HU  I  ii.".  II  I  |  •*  ‘ :  I  •  I 

■Tv)  DAT  A  I  t  in  ITVtIHl  I  M.III  null  I.  4  I  " s i •  > .  i  llv'VI  Ii  I.  Hll'l 

"'.0  DATA  IhAlilNAD  v  DAI  I  3 

60  DAT  A  H  > 2 

70  DAT  A  1 . 7.05  ,  .3 .  5 .  .  8  ’ .  2  . 5  •  .  5  >  I  .  5  •  .  I  ft  .  . ' . .  .  0  ", .  .  5 

80  DATA  .4rl.5r  .625,  2.5.  .835.  VS .  v, 

90  DATA  14.  '.4.44 

1 00  DA  T  A  0 .0.1  .  .  OA  .2.1,7;'-  S:i ,  4 .  .>  .44.  5  -  4. 3 .6.5.1  7 

1  10  DATA  7 .6. 0.3 ,8.  A.  V.  V.  7.  >'7.  I  0.8. 64.  I  I  I  7.10  -  ,7.  13.11.  '4 

130  DATA  -.895 
140  DATA  7 

145  DATA  -3. 1  .  A.  -2,  I  .  185.  1 .  .  71  5, 0 .  .  1 , 0 .  .  755.  I .  .  A6 2 .  -.85 
150  DATA  7 

155  DATA  3.1.5.  2.1.085,  1  .  .015,0,0.0,  •  .  355 » -•  1  .  .  76  ,  -2 ,  .95 
160  DATA  7 

165  DATA  -3, 1.4, -2, .985,  1 , .515,0,  . 1 .0, -.455,-1 .-.8A.-2.-l .05 
170  DATA  1 
180  DATA  8 

1  85  OA T A  -3 » 1 .6,-2. 1 , 29 » •  ■  1 , .92,0, .55, 1 , .01 , 1 ,  . 25 

190  DATA  0,  .65,-1, . 83 

200  DATA  8 

205  DATA  -3,1.5, -2 ,1.19, •  1 , .82,0, .45, 1 ,- .09, 1 ,-.3S 
210  DATA  0,-. 75,-1, -.93 
220  DATA  8 

225  DATA  -3, 1 .4,-2, 1 .09,-1 , . 72,0, .35, 1 . -. 19, 1 ,-. 45 

230  DATA  Of- .85, -1 , -1 .03 

240  DATA  2 

250  DATA  0 

260  DATA  9 

265  DATA  -3, 1 .59, -2, 1 .29, -1 , . 98,0, .68, 1 , . 37,2,-.02 
270  DATA  2, -.57,1, -.79,0, -.93 
280  DATA  9 

285  DATA  -3, 1 . 49 , -2, 1 . 19, -1 , . 88 ,0, . 58, 1 , . 27, 2,  .  12 
290  DATA  2, -.67,1,- .89,0,-1.03 


Statement 

Number 


10 

20 


30 

32 

34 

36 


40 

50 


60 


70 

80 

90 


100 

110 


Rl,  R2 ,  N9 ,  R8 ,  R9 
Al,  A2,  FI,  F2 ,  N3 

Statements  10  and  20  contain  the  six  moire  device 
parameters  and  other  values  used  to  establish 
limits  for  subsequent  computations  and  procedures. 
Note  that  N3  has  a  value  of  3,  indicating  that 
profiles  are  to  be  obtained  along  three  sections 
of  each  frame. 

P$(l)  to  P$ ( 9 ) 

P$(ll) 

PS (12) 

PS  (13) 

Statements  32  through  36  are  descriptions  of  the 
section  lines  which  were  used  for  each  frame. 

S$,  D$ ,  P$(0) 

PS  (10)  ,  Y8 ,  Y9 ,  NO 

Statements  40  and  50  contain  descriptions  of  test 
target  conditions  and  a  specification  of  the 
number  of  frames  which  will  be  analyzed  for  the 
particular  test. 

PI,  Ml 

The  number  of  points  and  the  order  of  fit  desired 
for  the  film  magnification  correction  polynomial. 

P  (1)  ,  Q  (1)  ,  P  ( 2 )  ,  R  ( 2 )  ,  .  .  . 

P  (6 )  ,  Q ( 6 )  ,  P (7)  ,  Q (7) ,  .  .  . 

Data  for  film  magnification  determination 

PI,  Ml,  U(0) 

The  number  of  points,  the  order  of  fit,  and  the 
distance  to  the  first  frame  of  the  event.  These 
values  are  used  with  the  film  speed  data  to 
determine  the  coefficients  of  a  polynomial  which 
describes  the  framing  rate  of  the  camera. 

R (1 )  ,  U(l) ,  R (2) ,  U ( 2) ,  . . . 

R(8) ,  U (8) ,  R(9) ,  U ( 9 ) , . . . 

Data  for  film  speed  determinations 
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Statement 

Number 


130 


140 


145 


150,160 


155,165 


170,250 


Variable  Name/Description 


X0 

Distance  from  origin  of  moire  device  coordinate 
system  to  point  of  interest  on  target 

PI 

Number  of  data  pairs  for  first  section  in  pre¬ 
event  frame 

U<1) ,  R(l) ,  U(2) ,  R (2 ) , . . . 

Fringe  number,  measured  value  pairs  for  first 
section  in  pre-event  frame 

PI 

Number  of  data  pairs  for  remaining  sections  in 
pre-event  frames 

U  (1)  ,  R  ( 1)  ,  U(2),  R(2)  ,  .  .  . 

Fringe  number,  measured  value  pairs  for  section 
in  pre-event  frame 

N2 

Frame  numbers 


180,200,220 
250,260,280  PI 

Number  of  data  pairs  for  each  section  in  the 
various  post  event  frames.  Note  that  a  value  of 
zero  was  entered  in  statement  250  indicating  data 
were  not  available  for  this  section. 

All  other 
statement 

numbers  Q(l)  ,  P(l),  Q(2),  P(2),... 

Data  pairs  for  sections  of  all  post  event  frames. 


The  next  three  pages  are  printout  that  was  received  at  the 
the  terminal  during  the  run  using  the  sample  data.  Only  results 
of  the  analysis  of  certain  pre-event  data  are  printed  at  the 
terminal.  These  data  are  printed  to  permit  the  user  to  select 
best  fit  data  for  use  during  analysis  of  subsequent  post  event 
frames.  In  each  of  the  pages  the  user  was  asked  for  a  response 
after  results  of  the  first  curve  fit  were  printed.  As  shown. 
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the  response — YES — resulted  in  another  printout  and  question. 

The  negative  response  resulted  in  a  request  for  specification  of 
the  order  of  fit  to  be  used  with  the  data  that  were  presented 
at  the  beginning  of  the  printout.  A  similar  sequence  was 
followed  for  the  additional  sections  of  the  pre-event  frame. 

At  the  end  of  the  third  page,  the  number  of  frames  which  had 
been  processed  and  a  notification  that  the  job  had  been  com¬ 
pleted  are  printed  to  signify  that  work  on  the  run  has  been 
terminated. 
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SAMPLE  TERMINAL  PRINTOUT 


First  Section 


*RUN 

>  USED  AS  ARRAY  &  SCALAR  NAME  2 E r H. M ,N , W 


SECTION  1-INCH  ABOVE  IMPACT  POINT 


FRINGE 

MEASURED 

X » MM 

DISPLACEMENT  »MM 

NO  TIMES  IN 

NUMBER 

RAU 

CORR 

RAU 

CORR 

SUBR  COR 

-3 

101.88 

100.86 

-24.02 

-25.41 

1 

82.30 

79 . 80 

-14.58 

-17.05 

2 

-  i 

53 . 76 

50 . 98 

-8.80 

-9.65 

2 

0 

6.21 

6.84 

1.14 

1.26 

n 

0 

-26 .51 

23.47 

4.88 

4.32 

3 

-l 

68.54 

68.37 

0.35 

0.26 

1 

•  89 . 99 

-93.41 

-6 . 72 

-3.23 

1 

COEFFICIENTS  OF  A  2  ND  ORDER  POLYNOMIAL  FIT  TO 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A  ( () )  =  1  .65370E  F00 

Ad)  =  -1  .03775E-01 

A ( 2 )  *  - 1 . 67575E -03 

STB.  ERROR  OF  ESTIMATE  ~  8.38858E--01 

THE  PRE  IMPACT 

FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X  t  MM 

DISPL »MM 

DISPL  r MM 

ACT  CAI  C 

-3 

10. 086E+01 

-25 . 4073E+00 

•25 . 859 IE  FOO 

45. 1782E-02 

■2 

79 . 800E+00 

-  1 7 . 0539E  TOO 

17. 3286E+00 

27 . 4673E  -07 

-  1 

50. 983E+00 

-96.4530E  01 

79.9291E -01 

-16.5230E--01 

0 

68.370E  Ot 

12.5744E  01 

86 . 585 5 E  02 

39. 1503E-07 

0 

-  2 3.47 2EL 00 

43. 1690E  01 

31 .6630E  01 

11 . 5061 E -01 

1 

-68.372EE00 

25 . 5676E  02 

9.1  .5397E-  07 

-65.9720E- 07 

-  2 

-93 . 414E  TOO 

-37.3170E  01 

-  32 . 7503E-01 

43.37A4E-03 

IT  TO 

NEXT  HIGHER  ORDER  POLYNOMIAL 

?  ( YES/NO ) ? YE 

CJ 

COEF  F  ICIENTS  OF  A  3  RD  ORDER  POLYNOMIAL  t  i  l  TO  I  HE  PRE  IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A ( 0 )  =  1.73528ET00 

A  <  1 )  =  l .  22948E- 01 
AC?)  =  -1.70283E-03 
ACT)  =  2.58886E  06 


STD. 

ERROR  OF  ES 

TIMATE  =  7.07723E-01 

FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X » MM 

DISPL »MM 

DISPIrMM 

ACT  CALC 

-3 

10 , 086ET01 

•25 • 4073E+00 

-25 . 3306E+00 

-76. 7212E -03 

-2 

79.880E100 

--17.0539E+00 

•  1 7 , 63 1 7E LOO 

57.7805E-02 

-  1 

50 . 9B3E+00 

96.4530E  01 

-86. 161 IE  01 

• 10.2919E-01 

0 

68.370E- 01 

12.5744E-01 

01 .5915E-02 

44.1524E02 

0 

- 23.472EL00 

43.1690E  01 

36.4950E  01 

66.7404E  02 

-1 

-68 • 372E+00 

25 . 5676E-02 

13.5380E-01 

10.9813E  01 

..  n 

-93.414E+00 

-32.3170E  01 

-  37. 4904E- 01 

51.7338E  02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  (YES/NO)?NO 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE 
PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  ?3 
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SAMPLE  TERMINAL  PRINTOUT  -  Second  Section 


SECTION  THROUGH  IMPACT  POINT 


FRINGE 

MEASURED 

X » MM 

DISPLACEMENT, MM 

NO  TIMES  IN 

NUMBER 

RAW 

CORR 

RAU 

CORR 

SUBR  COR 

-3 

97.64 

95.55 

■29.77 

-32.61 

■;> 

76.79 

72.57 

-20.20 

-24 . 35 

3 

-1 

46.82 

42 . 82 

-10.91 

-12.12 

3 

0 

■  2 . 62 

-  2 . 40 

0.48 

0.44 

1 

0 

36.42 

-32 . 02 

6 . 70 

5.89 

3 

-1 

79.70 

-75 . 68 

7.16 

4.08 

3 

..  '■> 

101.73 

-99 . 00 

5 . 26 

2.47 

■j 

COEFFICIENTS  OF  A 

2  NB  ORDER 

POLYNOMIAL 

FIT  TO 

THE  PRE-IMPACT 

PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 


A  <  0 )  ==  r-G .  451840-02 
A(1.)  =  - 1 . 92726E  01 
A (2)  =  - 1 . 68402E- 03 


STD. 

ERROR  OF  ES 

IT  I MATE  =  9 , 99634E  -01 

FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X ,  MM 

D I SPL, MM 

D I SPL , MM 

ACT  -  CALC 

-3 

95 . 546E+00 

32 . 6095E  TOO 

-33 . 8796E ) 00 

12.7014E-  01 

....  n 

72 . 566E+00 

-24 . 3455E+00 

-22 . 9420E+00 

1 4 . 0352E • 01 

~ 1 

42.825E+00 

-12. 1 L184E+00 

-11. 4279F 1 00 

-*69 . 049  IE-02 

0 

-24.050E- 01 

44.2338E  02 

36. 924 IE  02 

73.0975F -03 

0 

-32 . 01 6E+00 

58 . 8826E- 01 

43.5886E -01 

15 . 294 OE  01 

-1 

■  75 . 676E  TOO 

40.8426E  01 

40.51541::  0.1. 

76.7280E  02 

-90. 998EF00 

24.7134E-0J 

24.82791  01 

11. 4 44 IE -03 

FIT  TO  NEXT  HIGHER  ORDER  POI  YNOMIAl  ..?  ( YES/NO ) ? YES 


COEFFICIENTS  Of  A  3  RD  ORDER  POLYNOMIAL  FIT  TO  THE  PRE-IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A  <  0 )  - 1.351 81 E- 01 
Ail)  =  -2.20498E-01 
A<  2)  =  -1.66065E  03 
A(3)  =  3 • 73797E- 06 


STD.  ERROR  OF  ESTIMATE  =  7.54936E -01 


FRINGE 

MEASURED 

ACT  UAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X ,  MM 

D I SPL, MM 

D I SPL, MM 

ACT  -  CALC 

-3 

95 . 546E+00 

-  -32.6095E+00 

33. 1026ET00 

49 . 3179E-02 

-2 

72 . 566E+00 

-24 • 3455E+00 

■23 . 4523E+00 

-89 . 3219E-02 

-1 

42 . 825E+00 

-12 . 1 184ET00 

-12. 3299ET00 

21.1589E  -02 

0 

-  24 . 050E-01 

44 • 2338E-02 

38.5457E-02 

56. 881 7E -03 

0 

-32.016E+00 

58 . 8826E--01 

50.9945E- 01 

78 . B801E-02 

-1 

-  75.A74E  TOO 

40 . 8426E- 01 

54 . 2099E- 01 

-13.3A73E  01 

...  n> 

98 . 998ET00 

24.7134E-01 

17. 9169E-01 

67 . 9&55E-  02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  <YES/NO)?NO 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE 
PANEL  DISF'I  ACEMENT  FOR  POST  IMPACT  FRAMES  ***3 
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SAMPLE  TERMINAL  PRINTOUT  -  Third  Section 


SECTION  1-INCH  BELOW  IMPACT  POINT 


FRINGE 

MEASURED,  a. MM 

NUMBER 

RAW 

CORR 

-3 

93.10 

89.54 

-2 

70.98 

65 . 42 

-i 

39 . 57 

34.17 

0 

-11.75 

-10.57 

0 

-46. 44 

•40.57 

1 

91 . 15 

-81.51 

-"2 

-113.77  - 

104.38 

DISPLACEMENT. MM  NO  TIMES  IN 


RAW 

CORR 

SUBR  COR 

-35.92 

-40.75 

2 

-25.06 

-27.56 

9 

-  13.11 

14.74 

3 

2.16 

1.94 

2 

9.17 

7.46 

4 

19.61 

9.13 

p 

24 . 48 

7.97 

3 

COEFFICIENTS  OF  A  2  ND -ORDER  POLYNOMIAL  FIT  TO  THE  PRE-IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 


A ( 0 )  =  1.56133E+00 

A ( 1 )  =  -2.810B6E-01 

A( 2 )  =  1.803 

57E-03 

STD . 

ERROR  OF  ES 

TIMATE  =  6.90966E  01 

FRINGE 

MEASURED 

ACTUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X .  MM 

DISPL. MM 

DISPL .MM 

ACT  -  CALC 

-3 

89 . 540E+00 

-40 . 7533E  LOO 

41 . 1098EL00 

43.6462E-02 

..  n 

65 . 422E+00 

-27 . 5633E+00 

27.66 99EI00 

1. 0.661 3E -02 

-1 

34- 169E+00 

14. 7423EI 00 

13. 2713E+00 

•14. 7100E-01 

0 

- 10.566E+00 

19. 432 IE  01 

1 2 . 072 IE  01 

73.5997E  02 

0 

■  40 . 567E  LOO 

74 . 6073E-01 

68 • 7330E-01 

58 . 7427E-02 

-  1 

-81 .507E+00 

91 . 2654E  01 

93 . 6728E- 01 

-24 . 0739E- 02 

....  p 

-10 . 438E+01 

79. 7279E-01 

01 . 277 3E  01 

- 15.4938E- 02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  ( YES/NO >?YES 


COEFFICIENTS  OF  A  3  ED -ORDER  POLYNOMIAL  FIT  TO  THE  PRE-IMPACT 
PANEL  DATA  (REFERENCED  TO  THE  MOIRE  DEVICE). 

A  ( 0 )  =  1  .71  lit:’,  Eto  0 
A ( 1 )  =  -2.99229E01 
A<2>  -  ' 1 . 74 33 IE -03 
A ( 3 )  -  2.52445E-06 


STD. 

ERROR  OF  ES 

TIMATE  =  5 . 28900E- 01 

FRINGE 

MEASURED 

AC IUAL 

CALCULATED 

DIFFERENCE 

NUMBER 

X » MM 

DISPL. MM 

DISPL » MM 

ACT  -  CALC 

-3 

89.540E  LOO 

-40 . 7533FI 00 

40.6694E+00 

-83 • 8928E -03 

„  •;> 

65 . 422EL00 

-27 • 5633E  f 00 

28 . 04251+00 

47.9263E  02 

1 

34. 169E+00 

—  14 ~ 7423F  LOO 

-  1  3 . 8706E 1 00 

•  87. 1632E-  02 

0 

-  10. 566EE00 

1  9. 4321 E  0.1 

12.5216E  01 

69,  .1  051.  E  -  02 

0 

■40 . 567E+00 

74.6073E-0) 

73.8944E  -01 

71 .2891E-03 

-  1 

-81 . 507E+00 

91 • 2654E-01 

97 . 2894E-01 

-60 . 2397E-02 

-2 

--10.438E+01 

79 . 7279E- 01 

76 • 5639E-01 

31.6407E  -02 

FIT  TO  NEXT  HIGHER  ORDER  POLYNOMIAL?  ( YES/NO >?N0 

DESIRED  ORDER  OF  FIT  FOR  USE  IN  EVALUATING  RELATIVE 
PANEL  DISPLACEMENT  FOR  POST  IMPACT  FRAMES  ?3 


1  2  JOB  COMPLETED -  GOING  TO  'END'  STATEMENT 

* 


The  next  four  pages  present  the  sample  printout  produced 
at  the  high-speed  printer.  Presentation  of  these  pages  completes 
instructions  for  the  use  of  the  moire  device  and  program  PANDIS. 


f 1  fa  -  A  BIRDSTP1KE  PROGRAM 


shot  nuketc* 


TEST  PROnLEM 


DATE  FjPED 


ocTOdr.'  4  i9ao 


transparency 


CUKVlP  LINES 


P"CJECT  iLt 
»C lEHT 
VELOCITY 


IMAGINARY  dALL 
C.uC  Lf-S 
3  KNOTS 


FRAMES  of  FILM  (ATTEP  IMPACT*  LwAlUATLC  FOR  THIS  SHOT 


FOLLO.ING  APE  THE  COEFFICIENTS  OF  Lf  AS  T -SOU AWE  S  POLYNOMIALS 
uF  THE  FORM: 

Y  :  AC  ♦  »1*X  ♦  ...  ♦ 

FIT  10  THE  EXPERIMENTAL 


-AlA.  FILM  MAGNIFICATION,  IN  Pl*L  INCHES,  Is  GIVEN  AS  A  FUNCTION 
OF  MEASURED  INCHtS  AND  FILM  Sr  E  r  0  »  IN  PICTURES  PER  SECOND,  IS 
oIVEN  AS  A  FUNCTION  OF  FPAKf.  NlJMbrn. 


COEFFICIENTS: 

A  C 

A  1 

A2 

STD  CR  OF  ES 

FILM  MACNIF  : 

-t.U2°7SE-Jl 

3.r IS 12E*un 

-F.'VMURL-Cl 

7.UA2E3E-02 

FILM  Sr  EE D  : 

7.  s6027E«OI 

I  .  2  A  S  5  A  C  ♦  0  0 

-0.7897  It -02 

6. 92820E *00 

FLS  STf  1  4d  LOCATED  °5.24  MM  TO  7 HE  °  1GHT  OF  CENTER  OF  NOIRE  DEVICE 

**  RESULTS  OF  CURVE  FITS  To  PPE-l“P&CT  PANEL  SECTIONS  ♦♦ 
(REFERENCED  TO  THE  MuIRl  DEVICE! 

sEC  COEFFICIENTS  OF  L  E  A  S  T -S  >.U  A  *  t  POLYNOMIALS  STP  ERROR 


NO 

A  ( 0  > 

A  ( 1  > 

A  (2  > 

A  1  3) 

A  1  4  1 

OF  EST 

1 

1. 7J5E»r.C 

-1  .2.9F.-01 

-1.703P-03 

2 .589E-06 

O.DOCf *00 

7.08E-D1 

c 

-1. 3S2E-C1 

-2 .2CiSr -Cl 

-1 .fcfelE-n3 

3.736E-C6 

G . CO  CE ♦ CD 

7.5SE-31 

-1. 712E«00 

-2.9V2E-01 

-1.743E-C3 

2.52<*e-r6 

0. CODE  *U0 

5.29E-01 

PROGRAM  RUN  At  161105  hRS  ON  1105*0. 


SAMPLE  PRINTOUT  -  Title  Page 
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SHOT  NUMBER 


TEST  PPORLEM 


PAGE  2  OF  4  PAGES 


FRAME  NUHeEP  --  PRC  -  .IMPACT 


X-C00R0INA1E 

,  CliPL* CEMENT  , 

X-CoGPDINATE , 

DISPLACEMENT , 

r  RINSE 

IN  PK, 

IN  PM, 

1  'i  MM, 

IN  MP, 

NUMBER 

3/P  TO  MOInt 

A /R  TO  “OISE 

»/P  TO 

»/ R  TO  PPE SHOT 

CCVICE 

DEV  ICE 

FUS  STA  1 AO 

POSITION 

SECTION 

1-INCH  AtOVt 

IMPACT  POINT 

-  » 

10J. 653 

-25.A07 

196.095 

-0.077 

79.88 j 

-17.CSA 

175.117 

0.578 

-  1 

sc.ve3 

-9 . 6  A  a 

1A6.220 

-1.029 

n 

6  •  1 7  7 

1 .257 

102.07a 

w  .  A  A  2 

-23.A72 

A  .  317 

71.765 

•0.66  7 

~\ 

-6c. 272 

0.256 

26.865 

-1.09b 

w  •y 

-93.A IA 

-3.23? 

1.82? 

u.517 

SECTION 

THRCLGh  IMPACT  POINT 

-  3 

96.6U6 

-32.610 

1 9C . 7b  3 

0  •  A  9  3 

-? 

72.666 

-2  A  .  ?A  6 

167.803 

-C.893 

-1 

A2.t25 

-12.118 

1 3  6  .  r6  2 

0.212 

-2.A05 

r .  a  a  2 

92.P32 

0  •  C  r  7 

-32.016 

5.63? 

63.221 

u  •  789 

-1 

- 75  *  t  76 

A  .C  bA 

19.561 

-1.337 

-98.993 

2  .  A  7 1 

-3.761 

D.beo 

SECTION 

1 -Inch  PtLG* 

IMPACT  POINT 

-  X 

89. SAC 

-AO. 75? 

16 A . 77  7 

-O.OP.A 

*  7 

6S.A22 

-27.56? 

lb  0. 659 

0  •  A  7  9 

-1 

3  <(  .  1  6  9 

-  1 A  .  7  A  ? 

129. AOE 

-C.B72 

r* 

- 1 1!  .  666 

1.9A3 

oA  .67  1 

G.691 

L 

-AG. 667 

7.A61 

5  A  .  6  7  0 

0.071 

-J 

-el.5C7 

9.127 

13.73C 

-0.602 

-  2 

-10a . 3?A 

7.97? 

-9. 1A7 

w  .  3 1  6 

SAMPLE  PRINTOUT  Cont'd  -  Pre-event  Frame 


SHOT  NUMBER  --  TEST  PROBLEM  PAGE  3  OF  4  PAGES 

FRA  H£  NUMBER  --  1  I  C.347Q  MSEC  AFTEP  IMPACT) 


X-CCORDINATE , 

DISPLACEMENT, 

X -COORDINATE  , 

DISPLACEMENT , 

FRINGE 

IN  MM, 

IN  MM, 

IN  MM, 

IN  MM, 

NUMBER 

W/F  TO  MOIRE 

W/ft 

TO  MOIRE 

«/P  TO 

b/R  TO  PRF SHOT 

DEVICE 

DEVICE 

PUS  STA  ImD 

POSITION 

SECTION 

1-INCH  AbOVL  I 

MPACT 

POINT 

.  7 

ICC.tSo 

-2f .“07 

196 . 095 

-0.077 

_  -> 

86.57b 

-1C  .?1M 

161 .615 

9.779 

-1 

66.663 

-0.661 

lb?.  120 

12.469 

3 

MS. 006 

0 .?77 

1M  0.  ?4  3 

15.289 

1 

10.S21 

21 .T  37 

1US.756 

2L>.781 

1 

-12  .  ISM 

19. “b? 

o  3 . OB  3 

16.509 

'■ 

-SM .031 

13.900 

O.FD6 

10.953 

-1 

-79. *55 

7.MM1 

15.26? 

t.ces 

SECTION 

THROUGH  IMPACT 

POINT 

-3 

95.546 

-32.610 

190.763 

„.45i 

-  p 

PC. 216 

-16 .709 

175.465 

5.671 

- 1 

59. L  1  1 

-  6  .  c  6  7 

1 SM . ?M  7 

11.17b 

37.303 

6.660 

132. rM0 

17.33c 

l 

<  .  G  1  S 

<0.317 

■»  7 . 2S 1 

20.9C3 

l 

- 1  v .  b  ro 

2 1 .  m  a  7 

75.627 

1  7 . 8  0  ;> 

-61 .u9C 

1  7 . 7  2  2 

3M.  147 

11.437 

- 1 

-E  5. <,91 

1 3.020 

10.146 

6.720 

j  E  C  T  i  t  N 

1-INCH  cELO.  I 

Mp  ACT 

PC  INI 

_  7 

8</.b4j 

-MO. 7b7 

1  o  4 . 7  7  7 

-  U  •  0  5  4 

72.573 

165.170 

7 . 6  0  1 

-1 

51.372 

-9. '27 

146.609 

11.61b 

1 

29.257 

5  .  ?Bl 

12  4 .49  4 

1 7.27o 

1 

-6.517 

1 r  .  r  6  1 

b  6  »  7  2  C 

15.207 

1 

-27.051 

2  3 . 7  2  P 

oP.  166 

1 0.271 

-67. SIS 

2  1  .  r  2  <* 

27.722 

11.757 

-1 

-  °  •  !j  1  o 

1° . 7oM 

5.227 

9.107 

SAMPLE  PRINTOUT  Cont'd  -  Frame  #1 


SHOT  NUMBER 


TEST  PROBLEM 


PAGE  9  OF  9  PAGES 


FRAME  NUMBER  --  2  <  C.6938  MSEC  AFTER  IMPACT) 


X-COORDINATE  , 

DISPLACFMFNT  , 

X-COQPDINATE, 

DISPLACEMENT, 

FRINGE 

IN  MM, 

IN  “K, 

IN  MM, 

IN  MM, 

NUMBER 

W/P  TC  MOIfiL 

W/R  TO  MOTRF 

«/P  TO 

«/P  TO  PRfSHOT 

OEkICF 

DEVICE 

FUS  STA  mo 

POSITION 

CATA  N  f  T 

AVAILABLE  F u° 

SECTION  1-TnCH 

ABOVE  I "P  AC T  POINT 

SECTION 

THRCuch  impact 

PCIM 

-  • 

100.37? 

-it .rtf 

195.609 

9.151 

_  ~> 

36.578 

- 1  r .  ?  i  «• 

161. PIS 

19.033 

-1 

71 .CR9 

2.917 

166.336 

25.281 

r* 

55.796 

19 .589 

151 .033 

31.596 

l 

35.523 

25 .93T 

139.760 

37.197 

n 

c 

1 7.SC9 

it  .c71 

112.796 

91.C56 

*. 

-25.296 

37.927 

69.949 

33.609 

l 

-S2.o99 

2°. 96  5 

92.393 

2  3 . 65  7 

r 

-72.SP7 

29.525 

22.650 

lc.635 

SECTION 

1-INCH  PtLO.  IMPACT  POINT 

•  ! 

59.5*1 

-33.375 

190.218 

1  C  .  3  2  2 

-  ? 

8  C  •  2  1  6 

-16.7C5 

175.965 

16.922 

-1 

fc  3  .  fc  1  0 

-'*.251 

159.097 

23.997 

w 

9  7  .  c  1  9 

5.751 

192.556 

29.352 

1 

31.192 

I2.lt! 

126.928 

35.952 

IG.lD 

36 . c  7  1 

1UC .286 

91.963 

? 

-31 .51o 

99.Cfcr' 

63.321 

39 ,C°9 

1 

-5c. 518 

3  3 . 6  c  3 

36 .719 

29 »  36  C 

7 

-71.771 

2  ®  •  1  e  * 

16.966 

18.377 

SAMPLE  PRINTOUT  Concluded  -  Frame  #2 
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DATE! 
ILMED 

✓  . 

•/ 


